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ABSTRACT 


This study is to understand the effects of nonlinear 
behavior of structures using a finite element method. 

The nonlinear behavior equations are derived from 
equations of motion and constitutive equations. The basic 
theory is the principle of virtual work. 

The thesis begins with a comprehensive formulation of 
continuum-based finite element theory. The theoretical 


portion concludes with the details of both spatial and 


temporal discretization, including a discussion OF 
nonlinearity. 
In particular, large displacement problems and visco- 


elastic problems remain a challenging engineering problems 
today. The viscoelastic problems depend on the relaxation 


function which is the source of material nonlinearities. 
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PeEANALYots OF STRUCTURES 

Analysis of structures subjected to large displacements 
requires the inclusion of nonlinear terms in the strain 
displacement relations. In the formulation of the problem, 
either total Lagrangian or Eulerian approach is adopted. In 
order to assure structural integrity for loads beyond linear 
range, nonlinear analysis is very much mandated. 

Another class of problems is the viscoelastic behavior 
of structures - problems where stresses in the components 
vary with time at a constant deformation state or strains 
vary with time at constant stress states. 


These problems as applied to simple one-dimensional 


structures will be addressed in this research work. In what 
follows is a brief survey of literature in this area. At 
the end of this chapter, overview of the thesis is 
described. 


B. LITERATURE REVIEW 

A method for nonlinear static and dynamic analysis using 
cable and beam elements is given by Bergan et al {1]. 
Formulation allows for material nonlinearities and very 
large deformations. Gadala and Oravas {2} present a 


comparison of finite element formulations based on the 


Lagrangian, updated Lagrangian and Eulerian formulations for 
nonlinear mechanics problems. 

Henriksen [3] applies a single integral constitutive law 
to describe viscoelastic behavior in his formulation. 
Ghoneim et al [4] treat the total strain rate as sum of 
viscoelastic and viscoplastic contributions in implementing 
the constitutive equations in their formulation, Pinsky and 
Kim [5] use variational approach for studying the nonlinear 
viscoelastic shell behavior by representing the layered 
shell as a multi-director field. The formulation includes, 
within individual layers, the effects of transverse shear 
and transverse normal strain to arbitrary orders in the 
layer thickness coordinate [6, 7]. 

Saigel and Yang present a quadrilateral shell element 
and describe implementation of elastic-viscoplatic 
constitutive law for large deformation. A large collection 
of papers on the application of ADINA may be found in [8]. 

Response of nonlinear mechanical systems to random 
excitation is given by To [9]. Thermomechanical response of 
Uniaxial bars with thermoviscoplastic constitutive laws is 
discussed by Allen [10]. Large deformations of the 
contacting bodies on the contact surfaces and solution 
methods are given consideration in fll, ile Large 
deformations of composite beams subjected to axial dynamic 
load is reported in [13]. 


Based on Eulerian-Lagrangian formulation, frictional 


contact analysis is discussed by Haber et al [14]. 

Halleux and Casadei [6] present transient large strain 
analysis of structures based on biquadratic elements. 
Polar-decomposition is used in arriving at explicitly 
derived tangent stiffness matrices for a frame element by 
oem@oh and Atiuri [15]. A method based on co-rotational 
coordinate system, where the local cartesian coordinate 
system attached to each element translate and rotate is 
described by Mattiasson [16]. 

Inelastic problems based on incremental virtual work 
equations is given by Nagtegaal [17]. A constant strain 
rate is assumed during an increment in integrating the 
constitutive equations. 

Planar curved beams undergoing large rotations is 
analyzed using Helinger-Reissner mixed variational principle 
by Noor [18]. Steele discusses several nonlinearities 
including buckling, creep, and nonlinear elastic problems in 
f19]. The use of Lagrangian formulation using an elasto- 
plastic model is presented by Voyiadjis and Navaee [20]. 
Hirashima et al [21] present an approach for the snap- 
through analysis of curved shell panels and determination of 
initial stress matrix. 

Importance of viscoelastic relaxation effects on a 
pellet are investigated by Basombrio [22], where it is 
observed that the duration of the ramp input influences the 


response considerably. Viscoplastic constitutive models are 


obtained explicitly, without resorting to inversion, “to be 
directly implemented ina finite element formulation [23]. 
An intrinsic time scale that can vary with time and position 
is used, together with means to control the numerical 
integration in the approximation of constitutive model by 
Chambers and Becker [24]. 

Chandra et al (25, 26] discuss applications of visco- 

plastic models in the metal-forming. Elasto-plastic models 
are used in analyzing unilateral contact and friction 
problems by Cheng et al [12]. Chung et al [27] suggest 
methods for varying parameters in guiding solutions for the 
Simulation of plastic forming processes. Liu et al zea 
propose formulation for path-dependent materials based on 
arbitrary Lagrangian Euler method and include an elastic- 
plastic wave propagation problem. 
Rheinboldt and Riks describe several algorithms based on 
continuation method and mesh refinements for the solution of 
nonlinear finite element equations. Ryu and Arora [29] 
propose an approach to solve nonlinear finite element 
equations based on substructures concept. 

Several papers pertaining to visco-elastic, visco- 
plastic behavior and their implementation is reported in 
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C. THESIS OVERVIEW 

The objective of this thesis is to review some of the 
literature available in the analysis of large deformation of 
structures, to develop the theory co include large 
displacements and viscoelastic behavior of one dimensional 
structures, and in the process, to understand the physics of 
the effects of nonlinearities. Chapter II presents 
theoretical formulation based on principle of virtual work 
and incremental load method. Descretization based on a 
Simplex element is given. 

Chapter III gives the algorithms used in the solution 
procedure. Full Newton-Raphson and Modified Newton-Raphson 
methods are described. 

Chapter IV describes some test cases and numerical 
results obtained. 

Finally, Chapter V provides conclusions and some remarks 


for future work in this area. 


val 


Il. THEORETICAL FORMULATION 


Ape OV Vora 

In this chapter we consider a nonlinear continuum 
equations for a rod element. This approach is based on the 
principle of virtual work [48]. The resulting equations are 
in convenient form for numerical solution via finite element 
methods [49]. This, of course, is the conceptual basis of 
the discretized-solid finite elements, introduced by 
J.T.Oden et al [50] in the elastic structures and nonlinear 


continua. 


B. NONLINEAR CONTINUUM EQUATIONS 

Discuss a rod element with the total Lagrangian Method. 
It may be proved that internal virtual work is equal to 
external virtual work. Namely a deformable system is in 
equilibrium if and only if the total external virtual work 
is equal to the total internal virtual work for every 
virtual dileolacement consistent with the constraints. 

1. Principle of Virtual Work 

In this section we prove that total internal virtual 

work is equal to total external virtual work. Consider the 
concept of minimum potential energy [51]. The condition 
is satisfied to equilibrium state and the change of 


potential with respect to the displacement 


2 


must remain stationary. And consider the load- 


displacement and stress-strain relationship of figures 


shown below. 


sisuke- : 


Figure 2.1 Load-Displacement Curve 
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O€ 
Figure 2.2 Stress-Strain Curve 


In General, the total potential energy is equal to the 


summation of strain energy and potential energy. 


ie, Tt = U + W eZ) 
where, tp : Total potential energy 
U : Strain energy (internal work) 
W : Potential energy of load system (external 
work) 


From Figure (2.1), 


|, | p du 
= = > JPy (2.29) 
where, Pi °: > 7° Dread 
ui : deflection under it load 


From Figure (2.2), 
. 
Ul = | | o 6G cde] 


vol o 


es | GO GaaG eo.) (229 


vol 


Substituting the equation (222) and (2.3) in equation Zee 


we get 


= | 3 “OC “Gadicvoie = espa (2.4) 
p vol 1 i 


Equation (2.4) can be represented by matrix form, the 


following is the matrix formula. 


‘ 


T 
mp = J,.5 {o} {€} d(vol) - { u} { PB} C25) 


V 
where, {o } a vector of stresses 
{€} a corresponding vector of strains 
{ui} : a vector of elemental nodal displacements 
{pi } a vector of elemental nodel forces 


Let us define [B], a matrix relating strains and 


displacements, 
ie, {S} = [| Bae fig (2.6) 
and define [D], a matrix relating stresses and strains, 


or (CO = eb aiannoe, G23) 
14 


Substituting equation (2.6) & (2.7) in equation (2.5), we 
get 
T t T 


Tp = J Cee Loe Bisatvol) fur; - {ui} {pi} (2.8) 


vol 
Now we consider the minimum potential energy, the condition 
is as follows: 


ditp 


i 
© 


d{ui } 





Zio) 


Taking derivatives with respect to fui} in equation (2.8) 


and applying to equation (2.9), we get 


: i 
oi } = J (3 | PS eaiGzeds) efi } (C27 0)) 


Vo 


From relating loads and displacements, we can get the 


elemental stiffness matrix [k]. 


a 
ie, rm) = | Poe eo Te ls atvol) (eZ lel} 
vol 


Summing all the elements of the entire structure, equation 


(2.10) becomes 


{P} = ({K] fu} C225 
where, {P} = 2 {Pi} 
Peale LK | 
Misi 2 411i} 


2. One Dimensional Finite Element 
Consider a shape function which is linearly function 
With geometric boundary condition with one dimensional 


i 


finite element. We define that the shape function, Vy is 


Vy (0) =a, + by & 


In order to derive the constants in the Snape function tem 
rod element with linear displacement field, use the 
geometric boundary condition. 


We may obtain the constants as follows: 














X2 
al =e _ 
1 Lo 
1 
b ee (2.13) 
1 Lo 
X1 
Ae ges 
2 Lo 
1 
b a 
2 Lo 


The relationship between displacement and shape function is 


given below 


M 
ua Ses aon UNeilg MeP, 2...) 
= W(x) u(t) + KYox) u(t) (2.15 
Another consideration is the viscoelastic materials. For 


viscoelastic materials, the stresses may be expressed in 


terms of relaxation functions as given by 


sity: | att - Es &) af 
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Oat ee er 1 CL) (cena) 


where, See eer ore re) : elastic response 


G(t)- relaxation function. 
Assume the element is in equilibrium and quasistatic state. 
The motion equation for a one dimensional finite element is 


fiven by J.T.ODEN [ref.6], 


| MY C1 Ay vu ) av =P (u ,t) 


where, Vo = AobLo : undeformed volume 
E : total generalized force at node N 
N 
On substituting for above shape functions and their 


derivatives, we may obtain 


mmediver cre 6) = VYob, [| Se - I] (1 + D. Bes se (OLN a) 
0 (625106) ) 


C. CONSISTENT NODAL FORCES 
Consider a rod element in two configurations: undeformed 
and deformed. 
Let the distributed forces per unit length in these 
configurations be q and q* respectively. The consistent 
load vector depends on both the magnitude of the load and 


the detormed state. 


af 


One dimensional rod element under distributed load awe 


showed as following figures. 


$004) q*c4, &) 





where, Lo: undeformed rod element length 
L : deformed rod element length 
qa(x,t) is distributed load per unit undeformed length 


qa*(y,t) is distributed load per unit deformed Tengen 


Figure 2.3 One-Dimensional Rod Element under Distributed 


Load. 
Consider the relationship between the deformed and 
undeformed coordinate. It may be written as 

¥UX> LY So te ee (251s 


where x: undeformed coordinate of a point 
y: deformed coordinate of a point 
u: displacement of a point within element 


The consistent nodal forces vector is given by 


P= J eEN dV =[ eRNdY 


In index notation, nodal forces may .be written as 


Pe [ pF udV = BEY AV. (2.19) 
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with Nee, 2 


where F : body forces per unit mass deformed 
Fo : body forces per unit mass undeformed 
p ; deformed density 
QO : undeformed density 


0 
melating q*(y,t) and q(x,t) is given by 
OR, te) = A GhsGzy is) (2.20) 
where, r Sey 7 Gx >: stretch 
Consider the undeformed body force per unit mass, then we 


may write 


- Wdy Ete 
°" dm ~ PA. - pA. 


Using equation (2.17) 





a= x + u (2.22) 
M  M 
dy 
ae + u C2eaas) 
dx M,x  M 


The undeformed body force per unit mass is then 


q* 


A 
F = (1 +u,x um) (2.24) 


0 o Ao 





The consistent nodal forces take the following form, 
[ We 
Py = J, at Vgx + J ati, ,u Wy dx (2.25) 


Substituting the shape function in equation (2.25), we may 


get 


ine 


Lo 
P, = (1 + buy) G* (a, + Dy cima (22265 


Now we consider some special cases of equation (2.26). 


One case is the uniformly distributed loads. 
we 


x, ¢ 


Lo 


Figure 2.4 Rod Element under Uniformly Distributed Loads 


The uniformly distributed loads are given 
qa*(y,t) = qo h(t) 
Substituting the loads in equation (2.26), we may get the 


consistent nodal forces. 


P, go h(t) i 
hea So ee ; (2.23 





where, bh (be 
dx 


Another case is the linearly distributed loads. 





Figure 2.5 Rod Element under Linearly Distributed Loads 
20 


Linearly distributed loads are given 
y 
a(x,t) = 4 (t) + [a (t) - a wi 
Lo 
Substituting the loads in equation (2.26) we may get the 


consistent nodal forces. 


P, (Lo - uy, + uy) 2 i q, 
{ ES } = Se oe (252.5) 
p 6 1 z q, 


ft SOLUTION METHODS 

In this section we derive the equilibrium equations of 
incremental method, Full Newton-Raphson, and Modified 
Newton-Raphson Method. 

It is based on total Lagrangian Method and _ use the 
Taylor Serier at a given equilibrium point. 

The element stiffness matrices and consistent load 
vectors are assembled to form global matrices in the usual 
manner. The prescribed displacements are implemented by 
replacing the diagonal element of the corresponding row by 
unity and the rest of the element of the rows and columns by 
zeros [49]. 

dl. neremental Method (INC 
Consider the incremental process, we can rewrite the 
equation of load-displacement relationship, 
(Ki ({Ahiti1 = {APh+1 + {APcH (2.29) 
where, i: the current configuration coordinates. 
{APc}i : the force unbalance 


eal 


{AP}i+1: an increment in the external load 


The method is described by the Figure 2.6 





where, ABC : The algorithm of equation (2.29) 


AB’ C’ : The algorithm of purely incremental 
process 


Figure 2.6. The Algorithm of Incremental process 


The approximate values of A, E, ©, is computed from p =rfu). 
To improve these values, take the iteration process at each 
load level with {AFP}i+1= 0, as is Figure 2.7. 

The iteration processes are two ways, one is few 
load levels with many corrective iterations in each, or 
another is many ioad levels wit few iterations in each. 
This 18 sometimes calied ‘incremental with one-step Newton- 
Kaphson metnod”™. 

ce. Fuli Newton-Raphsor Method (FNR) 

Consider the load versus displacement function FP = 


f(u). Take the Taylor series of P = f(u) 


fith+ Adz $cup + (4 AL + 


(2-50) 
aie 
Where, etree = twee Caneent stifiness at A 
du A 
A 
Eu) ene 
A A 
memes 2 /\ jo = ok 
A ' B 
We look for A such hia’ (ul ND = el SO 
A B 
SOTA ean C2ob) 
A \ B B 


The consequent steps are explained by the Figure 2.7. 





Figure 2.7. The Newton Raplison solution of the equation 
p= f(u) 


The above algorithms are below; 
update displacements : ut = Ua + ae 
use ul to get stiffness ki and load Pt 


find next increment ae from ki (A,) =o eel 


Eventually ua + A, + ee a ="Un to a close 
approximation. 

Then, the load can be increased from Pa to Pc, and the 
precess can be started again. This solution algorithm is 


€ 
as 


known as the Newton-Raphson Method. 


3. Modified Newton-Raphson Method (MNR) 


The Newton-Raphson method requires that tangent 


stiffness matrix [k] be generated and reduced for equation 


solving in every iteration cycle. This process is not 
economic. An alternative is modified Newton-Raphson 
iteration, also called “constant stiffness iteration". Here 


the same tangent stiffness matrix is used for several 


iterations cycles. The MNR is described by the Figure 2.8 





LA 
Figure 2.8 The MNR Iteration 


As compared with FNR, the MNR requires more iterations, 
but each is done more quickly because [K] is formed and 


reduced in only the first iteration. 


KE. ANALYTIC SOLUTIONS 


At first we define the concept of stresses [52] Eulerian 


stress, o , is defined as follows: 


24 


Nitya) 


o (y,t) =-—_~— (62 oO ) 
A(y) 
and is called True stress. Lagrangian stress, T, is defined 


as follows: 


Ce C2) 


and is called Engineering stress. 


Kirchloff stress, S, is defined as follows: 





N(y56) 
SR pe) Shs (62-34) 
(x,t) Ao (x) 
The relationship between strain r , and stretch is as 
follows: 
a 
r= %@()X - 1) (2.35) 
du 
ra 1 + (2.36) 
dx 
where, N : Applied force 


A : deformed area 


For the special case of small strain, both u, x and change 
in cross-sectional area are small. In this case, we can 
write 


Yr = u,= € 


C=) 


where, € is engineering strain. 
Thus in small deformation case, the Lagrangian stress T 


and engineering strain € can be used to replace the 


Zo 


and engineering strain € can be used to replace the 
Kirchhoff stress S and Lagrangian strain r 
1. Elastic Materials 
Consider a one-dimensional rod with varying 
cross- sectional area, Ao = Ao(x), in equilibrium under 


distributed external loading. 


ax 
<i 


ea 
a Oe ——s(T+ 4h ox) 
4 A. Az (+ de ax) 
(a) (b) 
where, f : external distributed load per unit length 
T : Lagrangian stress within rod 


Figure 2.9. One-Dimensional Rod under Axial Loading 


For static problems and elastic materials, both f andT are 


independent of time, 


Cig te 0 Ree 


> > ab wt 4 s 


oe * 
Se j 


f 
ome \ ae 


From figure 2.9(b), the static equilibrium of forces in the 
y-direction 2s 


(2.37 
5 F(x) 2 -TAg + fax +(T + (dT/dx) Ax) (Ao + (dAg/dx) ax) = 0 


Expanding in Taylor series, 


Zo 


——- ax = Q (2.38) 
For a continuous rod, tie ene Limit. 
++ CTA.) =o (eas 


Seber ieticion Of eqiations (2.34), (2.35), and (2.36) in the 


above (2d O) 
where, u(x) : displacement w.r.t. x 

U, x >: tirst derivative of u(x) 

U,xx : second derivative of u(x) 


Note that equation (2.40) gives the values of distributed 
force f (oer unit length) required to produce the 
displacement field u(x) of a one-dimensional rod once the 
boundary conditions, the stress-strain law, and varying area 
function Ao(x) have been specitied. Consider the rod with 


constant area Ao, equation (2.40) becomes 


f(x) =-Ag uxx [(1+u,x) **2(dS* /dr) + S°] (ane 


oome special cases of equation (2.41) is considered. One 


case is the specified displacements. 


a 





where, u(x) = "aomeaxecs x= 


Figure 2.10 Test Model of “Inverse Solution” 


Apply the boundary “contititon 1 *(*) — moma soc amre 


a 
u(x) = ajeaee ae x2 (2.42) 
Ziie 
io ay 1S specified, we may get the exact solution u(x) and 


the corresponding required distributed load f(x). 


a. Linear displacements field: 
u(x) = a,x 
u, x = a, 
u,xx = 0 


From equation (2.42), we may get 
f(xy cee 
Note that a linear displacement field cannot be produced 


With distributed. loads 


Heme oOlUadratic displacement field: 


eCxXe) weet oe = cl /2Lo ) x¢ 
u, x merce yl aay) Lio ) x 
he = I lee) 
- Consider the linear materials, Se = Er 


We may vet the distributed forces from equation (2.41) 
Pex) = = BAwi, xx ele etiex, eet (17) x 2 to 2 x) | (G2n 43) 
= ~- ConSider the nonlinear materials, 
Sse = Bir + B2r2 (quadratic form) 
From equation (2.41), we may get 
Mie omeemeo fot xx { (1 + ,%) [ti 4+ 2 E2 (u,x + % ue,x)] 
+ Ei (u,x + %u2,x) + E2 (u,x + %u2,x)] } (2.44) 


Another case 1s the prescribed loads. 





Figure 2.11 Test Model of Rod with Prescribed End Element. 


From equation (2.39), since no distributed loads exist, 


d 
---—- (TAo) = 0 with constant area. 
dx 
We may get 
Se eons Cant (245) 


The applied load is given 
Sas (2.46) 


- Linear material; 


$e = Kr = ECU ee ene 
From equation (2.46), we may get 
Poe i Aol Cte uz,x + % us, x) (2.47) 


- Nonlinear material; 
Se «se, 1 r Wee (quadratic fom) 
From equation (2.46) we may get 
P = Ao (1 + u,x) [ Ei (u,x + % 2 fe) + BSG + Bee 


(2.48) 


NO 


Quasilinean Viscoelastic umcbelitames 
The Viscoelastic materials has the more complicated 
problems than the elastic materials, since the external 


leads and displacements are time-dependent. 


UC UL.C*%) 






seo at aa ee (A) 
X=e ete 
j= $$ $$______________—_—-J 
Eo 
where, Ao, Lo >: Underformed Area and Length 
Uo (t) : time-varying lett-end displacement 
uLo(t) : time-varying right-end displacement 
Po (t) : time-varying leftt-end force 
FLo(t) : time-varying right-end force 


Figure 2.12 One-Dimensional Rod Subjected to Time-varying 
Mad Loads 


3U 


For quasilinear viscoelastic materials, 


t 
Ss=Si) — | ett 8 508 dE 
0 6 


(CZ tS) 


The Lagrangian stress is constant throughout the rod since 
the distributed loads is not applied. From equation (2.39), 
the stress is time variable only. 

Met) = a(t) =) A(t) s(t) (2.49) 
Emel thie internal forces produced with the rod is 

Mju Tho = Aog(t) (2.50) 
Gubstituting equation (2.15) in equation (2.49) yields 

Phi Gt) S (ot) 


SaGtt= é 
S*(r) — gS Dg (r) a Zo) 
o0 a 


it 





(1 + u,) 


iiewtieit hand Side of equation (2.51) is a time function 


Omiy . 


Theretore g(t) must also be a time function. 


du 
—_—— ee hoes (x,t) = h(t) 2m) 
dx 
Pimteerating W.r.t. xX at equation (2.52) 
Metro t y= = xri( Ly) UF UC (4) (2255 )) 
Where, C : integrating function to be determined by the 


boundary condition at one end. 
Assume that 


Checr eects) ) = 0 ( t.) 


therefore 
(x,t) = xn ty) eter aere (2.545) 
Mex (ke tS eee 


Substituting equation (2.54) in equation 2a ote veld 
3a -=) Ly 
Fre h@oslS” Pho + GOS - J {KE} SMS HE] 55 
= ee Ck) 


Once the G(t) and Se@(r) have been specified, the equation 


55) can be used with the boundary conditions to find the 


ES 


( 
exact solution in some cases. Consider the 3 element solid 
with linear material property. The prescribed displacements 
boundary conditions at both ends is described in Figure 


Cagle 


UL. () 
SS xX 
Ao 
Lo 
B Ger S =.) Gow tment tlo (tee aw) 
uy (hem tye uy, (t) (2. S95) 


Figure 2.13 Test Model of a Rod with Prescribed Time Varying 


Displacement at the End. 


The prescribed displacements are linear, 


U(x) tae bow) up, (t) (2e one 
UO, « Ce eee u(t) J lio SiGe) Zea) 


C2 
CO 


where, Uy, (t) =o ( t) 


Assume linear material property, Se = Er 
And for a 3 element rod, G(t)is defined below 


G(t) 


Gia cr.) e-~ Bt (2.59) 
where, a , 6; constants for material properties. 


Note that G(o) oat G(t ) = £ - a at GB = 0 


No relaxation occurs and the material exhibits elastic 
response. 


Specifying the prescribed displacement function u p(t). 
a. Ramp Function 
Tao hiree os ul o,t) = 0 


u(ho,t)=,u (ro = au (a. constant) 


From equation (2.57) and 2.58), we may get 
mex st) weamta/bo.) xt (2.60) 
ex) «= \a/lo) t 

Pao lolblnie eguiatbions  (2.9/ ).-.(2.59) in equation (2.55) 


We may get the reaction force, F. 











EAo at az t2 aa a Sane 
ee elds ee et poo he els 
Lo Lo 2 Lo (3 BLo 
aa a 
= Smee orem te = 2Bb + 2) }) (ees. Oly) 
ei hie 


Equation (2.61) describes the reaction force produced at the 
fixed end for the prescribed displacement 


u ft) = at 


ers 


b. Step Function (Relaxation Test) 
At: x = Lo, BYG@e¥s - Stier.) 
u(Lo,t) = tagi(t) 
where, l(t) :© Unit Step funectien 
From equation (2.57) and (2.58), we may get 


uLo 
U Cee) = x 1(t) (2. Gey 





Uu,x = Lice.) (2 S6ian) 


Substituting equations (2.57) =~ (2.53) and (2.62) 


equation (2.55), we may get the reaction force, F. 




















EAo uLo ULo 
F(t) = ( ele 1a as ) G(t) (lL + (————) 1 Ga 
Lo 2 Lo Lo 
(2.64) 
Consider the Kirchhoff stresses. 
From equation (2.34) 
uLo uULo 
So2 ean ae a TG Cr tia we ne G(t) (2.659 
Lo 2 Lo . 
UL oO ULo 
Whe resol (lest ) « Lagrangian strain due to 
: Lo 2 Lo 


a prescribed step displacement function ul I(t) ac tnewaee 


end. 
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III. PROGRAM IMPLEMENTATION 


A. OVERVIEW 
We are interested in static analysis. Borestatic 


analysis, load-displacement relationship is as following: 


eek) 1; 


where {P} : external loads vector 

{k] : constrained stiffness matrix 

{u} : nodal displacements vector 
If we know {P} or f{u}, then we can analyze the static 
problems by using elementary elastic theory. In general, 


the static analysis steps are shown as follows: 


Step 1 : Enter geometry of material properties, loads, 
et al. 
Step 2 : Calculate element stiffness matrices. 


Step 3 : Determine appropriate part of [{k] and {P}. 
Step 4 : If [k] or {P} is appropriate, then repeat step 2 
through step 3. 
etep Oo): If [k]} or {P} a not appropriate, then apply 
constraints. 
Step 6 : Determine {u}. 
Step 7 : Regenerate elements. 


Step 8 : Calculate element stresses. 
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B. SOLUTION PROCEDURES 
In nonlinear structural problems, to find the solutions 
we use the INC and FNR and MNR. In fact, it is hard to know 
in advance what increments of loads or iteration should be 
used to obtain a good approximation to the solution. The 
interactive method was used in determining the strategy for 
solution in program application. The strategy subroutines 
were developed for elastic or visco-elastic time-dependent 
structural problems. Here the initial state for incremental 
loading process was selected as a null displacement vector. 
Next we consider finite element method for elastic (FEMEL) 
and viscoelastic (FEMEV) programs. IN both programs, 
equilibrium was checked at each stage for solutions obtained 
from other INC or FNR and MNR. 
1. Algorithm Implementation 
Consider the algorithm of Newton-Raphson Method 
based on total Lagrangian Method. Full Newton-Raphson 
Algorithm is as follows: 
Step 1. Establish local coordinates X-Y for the element at 
hand. 
Step 2. Generate the element stiffness matrix [k] in local 
coordinates so that it operates on the local D.O.F. 
Step 3. Transform [k] to global coordinates so that it 
operates on global D.O.F. 


Step 4. Repeat step 1 through 3 until all elements have been 
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treated. 
Step 5. Compute displacements. 
Step 6. Compute internal forces, [k]{u}. 
Step 7. Solve [k] {Au} = {AP} 
where, U,,, = Uyt A uy 
Step 8. Test for convergence. 
If not satisfied, return to step 1. 
Where, convergence = ((A RTA R) Vy (A pt Ap))**0.5 < 


Beo01 


af 


oe 


Flow Chart o 












Input:geometry,mat.prop;loads 


PROBLEM TYPE: 
ELAST/VISCOELASTIC 
ELE. MATRICES 
STIFF/LOAD: 


GLOBASSEMBLY 


SOLN: FORWARD 
SOLUTION 


SOLN:BACKWARD 
SUBSTITUTION 


UPDATE DISPS 


y 
COMPUTE RESIDUAL 
FORCES 


NEW LOAD INC./UPDATE DISPS 
PRINT RESULTS 





ae 








DIAGNOSTICS 


We consider two numerical examples for elastic and 
Viscoelastic materials. 
mee ELASTIC MATERIAL 


1. Prescribed concentrated forces. 


— a 





Opecify the rod element. 





pee ee OU Lis 


ee Ore oe | 


ie ho? fic) 
Pee 06 fips /ine | 
Consider the stiffness matrix of rod element, [k]}. For 


the rod element, 


EA Ki 1 Rae 
fae (4.1) 
k21 k22 


Consider the boundary conditions, 
U1 oa) 
Fl ane) 


Sick 


thie 


CO 


P2 = 22,5000 ataagewn 


From loads - displacements relationship, 


EA 1 —] al - 
L L ] ust |p 
Solve tor the displacements, u2, 


Hence, we nay get 





Consider the same rod element before, and specified 
quadratic displacements, u(x), 
uUtxX)--= O02 xe == Oe 


The first derivative is 


wee = 0,02 0200 Zar. 


The second derivative is 


deu 
miata ONE 
dx2 


4V 


From equation ( 2.41 ) with the constant area, we may 


getthe distributed loads for any x. 





d2u du dse 
E( X}=-=- Ao -—~ [ (1 + - } 2 wenn = + Ge | 
dx2 dx ar 
a OK. 10-3 ox2 oeio x) to) (4.2) 
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B. VISCOELASTIC MATERIAL 


Consider the rod element below, 


_ Fil) Fa ¢(k) 


ow 


Pommdary Conditions; uf(QO,t) = uo (t) = O 


Wiicow t )eee= eo Ul Cee = “u Pie.) 


Lo Lo 
maeraneian Stress; T(x,t) = eg(t) = A(t) S(t) 
imcernal gForce; F(t) = T Av g(t) 


P 
Kirchhoff stress ) Sr tegen ene mae setters emanate (4335 


N(x tale nor*( x ) 


For 3 element solids, the relaxation function is 


i wee de a) Fea Ee OE (4.4) 
where, 6 4 and ( gare constants ot the material 
properties. 


Assumé unit step function for prescribed displacements. 
At that time, the boundary condition is 
u(o, t) = O 


wii bo, tb a Ue be = Up 1 ( ta 
Oo oO 
where, 1(t): unit step function 
We may get the displacements function with unit step 
fone tLon . 
x 
u(x, t) = = ubo L(t) (4.5) 
Oe 
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du UL o 
— = Gc (4.6) 
dx Lo 





And the Kirchhoff stress is specified 





S =e Soe 
Substituting equations (4.4) - (4.7) am equation Gia. 
EK Ao hie: ULo 
F(t) = —-———- Ubo 1 + G(t) 1 +—— 1(t) (4.8) 
Lo ZLo Lo 


Hence, we may get the Kirchhoff stress, S, 


os rhs GEE (4.9) 
LES UL o 

where, Yiour= ry dees 
Lo 2 Lo 


Lagrangian strain 








due to the prescribed step displacement function 
u i(t) at the rod end. 
Equation (4.9) describes the relaxation cases. From 


equation (4.8), we may get the internal forces. 
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V. CONCLUSION 


This study is directed towards understanding the 
nonlinear behavior of structures. 

The formulation is based on the principle of virtual 
work. We consider the one dimensional nonlinear continuum- 
based rod element that includes large displacements and 
viscoelastic effects. To get the solution we use the 
incremental process with full Newton-Raphson and modified 
Newton-Raphson Method. Numerical solutions agree well with 
exact solutions. 

Further studies may be extended to two dimensional 
structures, such as laminated composite plates. 
Consideration of the viscoelastic effects of composite 
panels at high temperatures is another possible area of 
study. 

An experimental study to verify the results of the 


present study is also recommended. 
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